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CHARACTERIZING PROPERTIES OF A GEOLOGICAL 
FORMATION BY COUPLED ACOUSTIC AND 
ELECTROMAGNETIC MEASUREMENTS 

Background of Invention 

Field of the Invention 

[0001] The invention relates generally to characterizing properties of a geological 
formation saturated with a liquid. 

Background Art 

[0002] A logging operation in a borehole consists in physical measurements 
allowing to deteraiine properties of a formation surrounding the borehole. The 
measurements allow for example to detect an oil reservoir within the formation. 

[0003] FIG. 1 schematically illustrates an example of a logging operation from 
prior art. A logging tool 100 is located in a borehole 101 penetrating a formation 
102. The logging tool 100 comprises a sensor 109, e.g. an acoustic transducer, an 
electrode etc. The sensor 109 allows to measure logging data. 

[0004] The logging tool 100 is lowered in the borehole 101 by a cable 104 and 
slowly raised by a surface equipment 105 over a sheave wheel 106 while the 
logging data is recorded. A depth of the logging tool 100 is measured using a 
depth gauge 107 which measures a cable displacement. 

[0005] The logging tool may further comprise signal generating means 110, e.g. an 
acoustic transducer, an electrode, etc. The signal generating means 110 allow to 
excite the formation 102 with an exciting signal. A response signal generated by 
the exciting signal, e.g. an echo, or an induced signal, may be measured at the 
sensor 109. 
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[0006] The logging data acquired may be analysed either in situ near the logging 

tool 100, or analysed by a data processor 108, or stored, for later analysis. 

[0007] Typically, the formation 102 may be a porous medium comprising a soUd 
portion and an electrolyte. Electrokinetic phenomena may occur in such a 
structure, 

[0008] FIG. 2 schematically illustrates an electrochemical solid/electrolyte 
interface. An electric double layer 206 may be formed at an interface 201 between 
an electrochemical solid 202 and an electrolyte 203. Charged particles, e.g. anions 
204, may be located at a boundary of the electrochemical solid 202. Consequently, 
oppositely charged particles, e.g. cations 205, may tend to agglutinate close to the 
boimdary of the electrochemical solid 202, thus forming the electric double layer 
206. The electric double layer 206 comprise a first layer 207 arid a second layer 
208. The first layer 207 is composed of adsorbed cations that are attached to the 
interface 201. The second layer 208 known as diffused layer is composed of jfree 
cations. 

[0009] FIG. 3 schematically illustrates an example of a porous medixraL The 
porous medium comprises an electrolyte 303 and a matrix 302, i.e. a solid portion 
of the porous medium. The matrix is essentially constituted of solid particles 307. 
Each sohd particle 307 is surrounded by an electric double layer 306, due to an 
agglomerating of charged particles, e.g. cations 305 close to boundaries of the 
solid particle 307. 

[0010] A relative motion between the matrix 302 and the electrolyte 303 generates 
a motion of the free cations 305 of the diffused layer 306. As a consequence, in a 

case of a moving in a single preferred direction, the free cations 305 are moved 
following the single preferred direction, thus generating a current. 
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[0011] If an electromagnetic field is applied onto the porous medium, the charged 

particles of the diffused layers 306 may be moved upon the electromagnetic field. 
The solid particles 307 are hence also moved. 

[00121 When provided by an electromagnetic field, the moving of the solid 
particles 307 is referred to as an electromagneto-seismic effect in the present 
document 

[00131 Similarly, when provided by a seismic displacement, a displacing of the 
solid particles 307 or a flow of the electrolyte 303 is referred to as an seismo- 
electromagnetic effect in the present document 

[0014] Both phenomena are said to be electrokinetic phenomena. 

[0015] Various methods of prior art allowing to characterize a formation are based 
on the electrokinetic phenomenon. 

[0016] U.S. Pat, 2,814,017 to Schlumberger Well Surveying Corporation, 

published November 19, 1957, describes a first method for detecting permeable 
formations among a formation surrounding a borehole. The first method is based 
on the fact that if periodic pressure oscillations are passed through a permeable 
formation, a fluid filling interconnected pores of the permeable formation is forced 
to move in an oscillatory manner. The moving of the fluid is controlled by a 
viscosity and an inertia of the fluid. An effect of the viscosity on the moving of the 
fluid is a function of a permeability of the formation. For a plurality of formations 
of different permeabilities but containing fluids of a same viscosity and inertia, a 
difference in phase between an electric potential and a pressure hence varies with 
the permeabilities of the formations. 

[0017] In one embodiment, periodic pressure oscillations are produced in the 
formation. Indications are obtained of a phase displacement between the periodic 
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pressure oscillations and any periodically varying electro-filtration potential 
produced by the periodic pressure oscillation. 

[00181 ^ an other embodiment, a periodically varying electric field is established 
in the formation. Indications are obtained of a phase displacement between the 
electric field and any periodic pressure oscillation produced by the periodically 
varying electric field. 

[0019] In each embodiment, the phase displacement allows to detect the permeable 
formations and also to provide indications of their comparative permeabilities. 

[0020] U.S. Pat 6,225,806 Bl to Court Services Limited, published April 24, 

1997, describes a second method for characterizing a formation. A seismic source 
with two frequencies radiates radially a seismic signal within a borehole. A pair of 
electrodes above and below the seismic source allows to record a seismo-electric 
signal- The seismic source is lowered in the borehole and allows to provide a 
continuous logging of a formation surrounding the borehole. The recorded seismo- 
electric signal allows to provide a raw evaluating of the permeability of the 
formation. 

[0021] US Patent 5,417,104 to Gas Research Institute, published May 23, 1995 
describes a third method for determining a permeability of a sample formation. 

[0022] A pressure is applied to a fluid of the sample formation and an induced 
voltage signal iis measured using a pair of measurement electrodes near the 
application of said pressure. A differential fluid pressure in the sample formation 
is also measured between the electrodes. 

[0023] An electric potential is also applied to the sample formation. An applied 
voltage signal is measured using a pair of measurement electrodes near the 
application of the electric potential. An induced differential fluid pressure is 
measured between the electrodes. 
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[0024] An electrical conductivity of the rock is separately measured. 

[0025] The permeability is calculated as a simple ratio between applied quantities, 
i.e. the differential fluid pressure and the applied voltage signal, and induced 
quantities, i.e. the induced voltage signal and the induced differential fluid 
pressure. 

[0026] The third method may fail to provide a reliable measurement of the 
permeability of a formation surrounding a borehole. Indeed, a mudcake usually 
covers a wall of the borehole. 

[0027] US Patent 5,503,001, to Gas Research Institute published April 4, 1996 

describes a method allowing to remove effects due to a mudcake when applying 
the third method within a borehole. 

[0028] If a pressure is applied onto a wall of the borehole for a short duration, i.e. 
with an high frequency, a pressure drop occurs only in the mudcake. If a pressure 
is applied for a longer duration, i»e. with a low frequency, the pressure drop occurs 
both in the mudcake and in a formation surrounding the borehole. A characteristic 
frequency fs marks a maximal change with respect to frequency between the low 
frequencies and the high frequencies. 

[0029] The applied pressure diffuses into a fluid of a porous medium according to 
a diffusion equation. In time t the pressure drop penetrates a distance 1 according 
to the equation : 

[0030] F-Dt 

[0031] Where D is a diffusion constant, the diflfosion constant being characteristic 
of the porous medium. 

[0032] The method described above, when performed over a range of frequencies, 
allows to determine the characteristic frequency fg. As the diffusion constant of the 
mudcake is known, a thickness Icake of the mudcake may be evaluated. 
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[00331 US Patent 5,841,280, to Western Atlas International, published November 

24, 1998, describes a fourth method that allows to estimate a porosity of a 
formation surrounding a borehole, A transmitter imparts acoustic energy impulses 
into a liquid filling the wellbore* Hie acoustic energy in:5)ulses reach a wall of the 
wellbore and propagates in the formation. An acoustic response signal may be 
detected* An induced electrical voltage, i.e. a seismo-electric signal, is also 
recorded. 

[0034] The transmitter may generate various modes of acoustic waves within the 

formation, e.g. longitudinal waves, shear wave, Stoneley wave. 

[0035] A Stoneley induced electrical voltage associated with a Stoneley wave has a 
relatively large amplitude. A seismo-electric Stoneley signal is synthesized from 
the detected acoustic energy and from an initial value of the porosity. A difference 
between the synthesized seismo-electric Stoneley signal and the recorded Stoneley 
induced electrical voltage is calculated. The initial value of the porosity is adjusted 
according to the difference. The synthesizing of the seismo-electric Stoneley 
signal, the calculating of the difference and the adjusting of the porosity are 
repeated until the difference drops below a predetermined threshold, thus allowing 
to estimate the porosity. 

[0036] Such an iterative inversion method may be applied to an estimating of a 

conductivity or a porosity of a fluid of the formation. 

[0037] The fourth method, xmlifce the third method, takes into consideration a 
propagating of a seismic wave in a porous medium of the formation. 

[0038] A model of electrokinetic phenomena based on both MaxwelPs equations 
relating to electromagnetic phenomena and Biot's equations relating to the 
propagation of the seismic wave has been described by Steven Pride ("Governing 
equations for the coupled electromagnetic and acoustic of porous media", Phys. 
Rev., 1994). 
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[00391 In a case of aa acoustic wave having an angular frequency co, the frequency 

a)/27i: is within a sonic range or within an ultrasonic range. Electromagnetic signals 
Avithin the sonic range or within the ultrasonic range have relatively large 
wavelengths, which allows to keep only quasistatic terms in the Maxwell's 
equations. A coupling between an electric field and a magnetic field may be 
considered as low and xmsteady effects in the MaxwelPs equations are neglected, 

[0040] A resulting Pride's coupling equation may hence be written as: 

[00411 J-^0(S)E+Ii4^^^ (eq. 1) 

[0042] wherein: 

[0043] J represents an electric current density; 

[0044] a(a)) represents an overall electrical conductivity of the porous medium as a 

ftinction of the angular frequency o; 

[0045] E represents the electric field; 

[0046] L(o)) represents a frequency dependant electrokinetic coupling coefficient; 

[0047] P represents a fluid pressure; 

[0048] pj- represents a fluid density; and 

[0049] Us represents a displacement vector of the solid particles. 

[0050] The fourth method uses the coupling equation (eq. 1) for estimating the 
conductivity of the fluid of the formation. 

[0051] The third method fails to take into consideration a possible propagating of a 

seismic wave in a matrix of the porous medium. The coupling equation between a 
pressure and an electric field is written in such a streaming model as: 

[0052] J = ct((d)E - L(cd)VP 
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[0053] Evaluating a permeability of the porous medium is hence performed in the 
third method by a simple ratio calculation. 

[0054] Furthermore, a second coupling equation may be written ia the streaming 
model as : 



[0055] 



- io)W = L(a?)E VP 

77 



[0056] wherein W represents a relative fluid-solid motion; 

[0057] k((X)) represents a permeability; and 
[0058] r| represents a fluid viscosity. 

[0059] The model described by Steven Pride leads to a Pride's second coupling 

equation taking into consideration the possible propagating of the seismic wave in 
the porous medium: 



- icoW = L{co)E + [- VP + co^PfU^ ] 



[0060] -l^rv = i^yoj)!!. ^ p v ^ -r a/ fj^u^ j 2) 

Summary of Invention 

[0061] In a first aspect, the invention provides a method for characterizing a 
fomiation. The method comprises exciting the formation with an acoustic wave 
propagating into the formation. A seismo-electromagnetic signal produced by the 
acoustic wave within the formation is measured. The method further comprises 
exciting the formation with an electromagnetic exciting field. An electromagneto- 
seismic signal produced by the electromagnetic exciting field within the formation 
is measured. The measured seismo-electromagnetic signal and the measured 
electromagneto-seismic signal are analyzed to evaluate characterizing parameters 
of the formation. 
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[00621 In a first preferred embodiment, the acoustic wave and the electromagnetic 

exciting field are generated at a logging tool positioned within a borehole 
surrounded by the formation. 

[0063] In a second preferred embodiment, an acoustic response signal is measured. 
The acoustic response signal is produced by the acoustic exciting. Acoustic 

properties of the formation are estimated from the acoustic response signal. An 
electromagnetic response signal is measured, the electromagnetic response signal 
being produced by the electromagnetic exciting. Electromagnetic properties of the 
formation are estimated from the electromagnetic response signal. 

[0064] In a third preferred embodiment, the method further comprises selecting 
initial values of inversion parameters. A synthesis seismo-electromagnetic signal 
and a synthesis electromagneto-seismic signal are synthesized using the initial 
values of the inversion parameters. A first difference between the synthesis 
seismo-electromagnetic signal and the measured seismo-electromagnetic signal is 
calculated. A second difference between the synthesis electromagneto-seismic 
signal and the measured electromagneto-seismic signal is calculated. The values of 
the inversion parameters are adjusted according to the first difference and to the 
second difference. The method further comprises repeating the synthesizing using 
the adjusted values of the inversion parameters, the calculating of the first 
difference, the calculating of the second difference and the adjusting until the first 
difference and the second difference respectively drop below a first predetermined 
threshold and a second predetermined threshold. 

[0065] In a fourth preferred embodiment, the inversion parameters are an 
electrokinetic coupling coefficient and a mobility. The synthesizing is simplified 
by synthesizing only a synthesis seismo-electromagnetic slow longitudinal signal 
and a synthesis electromagneto-seismic slow longitudiaal signal from a mobility 
initial value and from an electrokinetic coupling coefficient initial value. 
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[0066] In a fifth preferred embodiment, the analyzing takes into consideration the 

propagating of the acoustic wave within the formation, 

[0067] In a sixth preferred embodiment the seismo-electromagnetic signal is a 

seismo-electric signal. 

[0068] In a seventh preferred embodiment, the seismo-electromagnetic signal is a 

seismo-magnetic signal. 

[0069] In a eighth preferred embodiment the electromagneto-seismic signal is a 
magneto-seismic signal. 

[0070] In a ninth preferred embodiment, the electromagneto-seismic signal is an 

electro-seismic signal. 

[0071] In a tenth preferred embodiment, the logging tool is displaced along the 

borehole so as to provide a continuous characterizing of the formation as a 
function of depth. 

[0072] In a second aspect, the invention provides a system for characterizing a 
formation surroimding a borehole. The system comprises a logging tool to be 
lowered into the borehole. An acoustic emitter located onto the logging tool allows 
to excite the formation with an acoustic wave propagating within the formation. 
An electromagnetic receiver allows to measure a seismo-electromagnetic signal 
produced by the acoustic wave within the formation. The system further comprises 
an electromagnetic emitter located onto the logging tool. The electromagnetic 
emitter allows to excite the formation with an electromagnetic exciting field. An 
acoustic receiver allows to measure an electromagneto-seismic signal produced by 
the electromagnetic exciting field within the formation. Processing means allow to 
analyze the measured seismo-electromagnetic signal and the measured 
electromagneto-seismic signal so as to evaluate characterizing parameters of the 
formation. 
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[00731 Iti an eleventh preferred embodiment, the electromagnetic receiver is an 

electric receiver allowing to measure a seismo-electric signal produced by the 
acoustic wave within the formation. 

[0074] In a twelfth preferred embodiment, the electromagnetic receiver is a 
magnetic receiver allowing to measure a seismo-magnetic signal produced by the 
acoustic wave within the formation. 

[0075] In a thirteenth preferred embodiment, the electromagnetic emitter is an 
electric emitter allowing excite the formation with an electric exciting field. 

[0076] In a fourteenth preferred embodiment, the electromagnetic emitter is a 

magnetic emitter allowing excite the formation with a magnetic exciting field. 

[0077] In a fifteenth preferred embodiment, the system further comprises at least 
one additional electromagnetic receiver and at least one additional acoustic 
receiver. 

[0078] Other aspects and advantages of the invention will be apparent from the 
following description and the appended claims. 

Brief Description of Drawings 

[0079] FIG. 1 schematically illustrates an example of a logging operation from 
prior art. 

[0080] FIG. 2 schematically illustrates an electrochemical solid/electrolyte 
interface. 

[0081] FIG. 3 schematically illustrates an example of a porous medium. 

[0082] FIG. 4 schematically illustrates an example of a system according to the 
present invention. 

[0083] FIG. 5 schematically illustrates an example of a system according to a 
preferred embodiment of the present invention. 
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[0084] FIG. 6 is a flowchart illustrating an example of a method according to an 
embodiment of the present invention. 

[0085] FIG. 7 illustrates an example of an algorithm for evaluating characterizing 
parameters of a formation surrovmding a borehole, according to an embodiment of 
the present invention. 

Detailed Description 

[0086] Acoustically exciting a formation generates a seismo-electromagnetic signal 

that comprises a seismo-electric signal and/or a seismo-magnetic signal. An 
electric field or a difference of electrical potentials may be measured, thus 
alloAving to measure the seismo-electric signal. Alternatively, a magnetic field is 
measured, thus allowing to measure the seismo-magnetic signal. Alternatively, 
both the electric field and the electromagnetic field may be measured. 

[0087] In the present description, the term "seismo-electromagnetic" may 

designate a seismo-electric signal produced by an acoustic signal or a seismo- 
magnetic signal produced by the acoustic signal. 

[0088] Similarly, it is possible to provide an electrical exciting of the formation : in 

this latter case, the electromagneto-seismic signal induced within the formation 
may be designated as "electro-seismic''. If the exciting is only magnetic, the 
electromagneto-seismic signal induced within the formation may be designated as 
*'magneto-seismic" . 

[0089] In the present description, the term ''electromagneto-seismic" designates 
either an electro-seismic signal produced by an electric excitation or a magneto- 
seismic signal produced by a magnetic excitation. 

[0090] Characterizing a formation surrounding a borehole, e.g. evaluating a 
permeability of the formation or a formation factor of the formation, allows to 
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assess a potential moving of a fluid in a porous medium. A quality and a quantity 
of water or oil reservoirs may subsequently be estimated. 

- [0091] Evaluating the permeability also allows to determine conditions of a further 
oil production. 

[0092] The formation parameter allows to evaluate a quantity of oil in the 

formation. 

[0093] The third method from prior art involves transient electrokinetic 
measurements that are rendered difficult by a mudcake of the borehole. 

[0094] The second method and the fourth method from prior art provide a 

measurement of a seismo-electric signal produced by a seisnciic signal. Such a 
measurement allows to evaluate a porosity of the formation, as performed in the 
fourth method. 

[0095] In the second method, the seismic signals are generated at multiple 

frequencies and amplitudes of the seismo-electric signals produced by the seismic 
signals are measxnred at each frequency, which allows to evaluate the permeability 
of the formation. 

[0096] There is a need for a method allowing to provide an accurate characterizing 
of a formation. 

[0097] FIG. 4 schematically illustrates an example of a system according to the 

present invention. Acoustic waves and electromagnetic fields are represented on 
FIG. 4 by arrows 408, 409, 410, 411, 412 and 413. 

[0098] The system allows to characterize a formation 401. The formation 401 
surroxmds a borehole 402. The system may comprise a logging tool 407 to be 
lowered into the borehole 402. Alternatively, the system may be located at a 
surface of the formation 401 . 
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[00991 An acoustic emitter 403 allows to excite the formation 401 with an acoustic 
wave 408 propagating within the formation 40L A seismo-electromagiietic signal 
409 produced by the acoustic wave 408 within the formation 401 is measured 
using an electromagnetic receiver 404. 

[00100] An electromagnetic exnitter 404 allows to excite the formation 401 with an 
electromagnetic exciting field 411. An electromagneto-seismic signal 413 
produced by the electromagnetic exciting field 411 within the formation 401 is 
measured using aa acoustic receiver 403. 

[00101] Processing means 414 allow to analyze the measured seismo- 
electromagnetic signal 409 and the measured electromagneto-seismic signal 413 
so as to evaluate characterizing parameters of the formation 401. The processing 
means 414 may be part of the logging tool 407 as represented in FIG. 4. 
Alternatively, the processing means are located at a distinct location, e.g. at 
surface. 

[00102] As both the seismo-electromagnetic signal 409 and the electromagneto- 
seismic signal 413 are measured, a method according to the present invention 
allows to provide a more accurate evaluating of a permeability of the formation 
40 1 than the second method and the fourth method from prior art. 

[00103] The seismo-electromagnetic signal 409 is produced by the acoustic wave 
408 that propagates within the formation 401, In the third method Scorn prior art, 
an applied pressure diffuses into a fluid of the formation. The permeability is 
evaluated by a simple ratio calculation based on a streaming model that fails to 
take into consideration a propagating of an acoustic wave in a porous medium. 
The method according to the present invention allows to provide a more accurate 
evaluating of the permeability of the formation 401. 

[00104] Furthermore, the method according to the present invention may provide a 
faster evaluating of the permeability or the other characterizing parameters than 
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the third method from prior art. In particular, the logging tool 407 may preferably 
be spaced from walls of the borehole 402 and may be easily displaced so as to 
allow a continuous characterizing of the formation as a function of deplii. 

[00105] Alternatively, the logging tool may contact the walls of the borehole. 

[00106] In the system of the present invention, the acoustic emitter 403 and the 
acoustic receiver 403 may be a single device, as represented in FIG. 4. However, 
the acoustic emitter 403 and the acoustic receiver 403 may also be distinct devices. 

[00107] One or more additional acoustic receiver(s) (not represented on FIG. 4) may 

be provided. The acoustic receiver and the additional acoustic receivers may form 
an array of acoustic receivers. 

[00108] Similarly, the electromagnetic emitter 404 and the electromagnetic receiver 
404 may be a single device, as represented in FIG. 4. However, the 
electromagnetic endtter 404 and the electromagnetic receiver 404 may be distinct 
devices. 

[00109] The system may further comprise one or more additional electromagnetic 
receiver(s) (not represented on FIG, 4). The electromagnetic receiver and the 
additional electromagnetic receivers may form an array of acoustic receivers, 

[00110] Preferably the acoustic emitter 403 and the electromagnetic emitter 404 are 
located onto the logging tool 407. The acoustic wave 408 and the electromagnetic 
exciting field 411 are hence generated at the logging tool 407. 

[00111] The measuring of the seismo-electromagnetic signal may be performed by 
measuring an electric field only, i.e» the seismo-electromagnetic signal is a seismo- 
electric signal. 

[00112] Alternatively, the measuring of the seismo-electromagnetic signal may be 
performed by measuring a magnetic field only, i.e. the seismo-electromagnetic 
signal is a seismo-magnetic signal. 
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[00113] Alternatively, both the electric field and the magnetic field are measured. 

[00114] Similarly, the exciting of the formation with the electromagnetic exciting 
field is only electric: an electro-seismic signal is subsequently measured. 

[00115] Alternatively, the exciting of the formation with the electromagnetic 
exciting field is only magnetic and a magneto-seismic signal is measured. 

[00116] Alternatively, the formation is excited with both an electric exciting field 
and a magnetic exciting field, either simultaneously or not. In this latter case, more 
signals may be measured, thus allowing to insure a reliable evaluating of the 
characterizing parameters. 

[00117] Preferably aa acoustic response signal 410 is measured ajEler the exciting of 
the formation 401 with the acoustic wave 408. The acoustic response signal 412 is 
produced by the acoustic exciting. The acoustic response signal 410 may be 
measured at the acoustic receiver 403, as represented in FIG. 4, or at any other 
receiver allowing to measure an acoustic signal. 

[00118] Similarly, an electromagnetic response signal 412 may be measured after 
the exciting of the formation with the electromagnetic exciting field 411. The 
electromagnetic response signal 412 is produced by the electromagnetic exciting. 
The electromagnetic response signal 412 may be measured at the electromagnetic 
receiver 404, as represented in FIG. 4, or at any other receiver allowing to measure 
an electromagnetic signal. 

[00119] FIG. 5 schematically illustrates an example of a system according to a 
preferred embodiment of the present invention. The system allows to characterize 
a formation 501 surroimding a borehole 502. The system coniprises a logging tool 
507 to be lowered into the borehole 502. 

[00120] An acoustic emitter 503 and an electromagnetic emitter 504 are located 
onto the logging tool 507. The acoustic emitter 503 and the electromagnetic 
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emitter 504 respectively allow to excite the formation 501 with an acoustic wave 
(not represented on FIG. 5) and an electromagnetic exciting field (not represented 
on FIG. 5). 

[00121] Preferably the acoustic exciting and the electromagnetic exciting are not 
performed simultaneously. 

[001221 The system further comprises an array of acoustic receivers 505 and an 
array of electromagnetic receivers 506 disposed onto the logging tool 507. 

[00123] Various modes of the acoustic wave may propagate within the formation, 
e.g. a fast longitudinal wave or a slow longitudinal wave, a shear wave, a Stoneley 
wave. The array of acoustic receivers 505 and the array of electromagnetic 
receivers 506 allow to properly measure signals providing from a desired mode. 

[00124] Furthermore, as a plurality of acoustic receivers 505 is disposed onto the 
logging tool 507, acoustic properties of the formation 501 may be easily 
determined at a further analyzing, e.g. by a simple ratio between a time of 
detection and a distance between two acoustic receivers 505. Similarly, 
electromagnetic properties of the formation 501 may be easily determined due to 
the array of electromagnetic receivers 506. 

[00125] By "electromagnetic receiver", we designate either an electric receiver or a 
magnetic receiver. The electric receiver and the magnetic receiver respectively 
allow to measure a seismo-electric signal and a seismo-magnetic signal produced 
by the acoustic wave within the formation. The electromagnetic receivers 506 may 
for example be electrodes, coils, magnetometers etc. depending on which field 
(electric or magnetic) is to be measured. 

[00126] Preferably the electromagnetic receivers 506 are relatively close to a wall of 
the borehole 502 so as to so as measure signals having a relatively high amplitude. 
The electromagnetic receivers 506 may also be relatively close to the 
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electromagnetic, emitter 504: typically a distance between each electromagnetic 
receiver 506 and the electromagnetic emitter 504 is relatively small compared to a 
wavelength of the acoustic wave. 

[00127] The acoustic receivers 505 may for example be made of a piezoelectric 
material. However, any other receiver allowing to measure a mechanical signal, 
e,g, a pressure field or particles velocities, may be used. Preferably the acoustic 
receivers 505 are relatively close to a wall of the borehole 502 so as to measure 
signals having a relatively high amplitude. 

[00128] If a vector quantity, e.g. a magnetic field, is measured, the electromagnetic 
receivers 506 may either allow measurements relative to a single axis, to two axes 
or to three axes. 

[00129] Similarly, if a vector quantity, e.g. a motion of particles, is measured, the 
acoustic receivers 505 may either allow measurements relative to a single axis, to 
two axes or to three axes. 

[00130] By "electromagnetic emitter'', we designate either an electric emitter or a 
magnetic emitter. The electric emitter allows to excite the formation with an 
electric exciting field and the magnetic emitter allows to excite the formation with 
a magnetic exciting field. The term "electromagnetic exciting field'' hence 
designate either the electric exciting field, the magnetic exciting field or a 
combination. Typically, a frequency range of the electromagnetic exciting field is 
such that a coupling between an electric field and a magnetic field in Maxwell's 
equations may be neglected: the electromagnetic exciting field is considered to be 
either the electric exciting field or the magnetic exciting field. 

[00131] The electromagnetic emitter 504 may for example be an electrode, a coil, or 
any other device allowing to create an electric field or a magnetic field. 
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[001321 The acoustic emitter 503 may be a transducer, or any device allowing to 
excite the formation 501 with an acoustic wave. The acoustic emitter 503 may be a 
monopole, a dipole, or a multipole. The dipole allows for example to excite the 
formation 501 with a shear acoustic wave. 

[001331 Preferably the electromagnetic emitter 504 and the acoustic emitter 503 
have a same frequency content. The acoustic wave and the electromagnetic 
exciting field have a same frequency so as to allow a further analysis of measured 
signals. Preferably the frequency is below a critical frequency of the formation 
501. Typically, the frequency is below 5 kHz. 

[00134] The acoustic exciting and the electromagnetic exciting allow to measure 
four signals (not represented on FIG. 5) : 

[001351 - an acoustic response signal, which is produced by the acoustic 
exciting- The acoustic response signal is measured at the array of acoustic 
receivers 505. 

[001361 - a seismo-electromagnetic signal produced by a propagating of the 
acoustic signal within the formation 501. The seismo-electromagnetic signal is 
measured at the array of electromagnetic receivers 506. 

[001371 - an electromagnetic response signal which is produced by the 
electromagnetic exciting. The electromagnetic response signal is measured at the 
array of electromagnetic receivers 506. 

[001381 - an electromagneto-seismic signal produced by the electromagnetic 
exciting field within the formation 501. The electromagneto-seismic signal is 
measured at the array of electromagnetic receivers 506. 

(001391 The measuring of the four signals allows to provide a further analysis. The 
further analysis allows to evaluate the acoustic properties, the electromagnetic 
properties and characterizing parameters of the formation 501, using for example 
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processing means 508. The processing means 508 may be a computer located at 
surface, a processor within the logging tool, or any other device allowing to 
evaluate the characterizing parameters. 

[00140] FIG. 6 is a flowchart illustrating an example of a method according to an 
embodiment of the present invention. The method allows to characterize a 
formation surrounding a borehole, 

[00141] The method comprises exciting a formation with an acoustic wave 
propagating into the formation (box 61). A seismo-electromagnetic signal 
produced by the acoustic wave within the formation is subsequently measured 
(box 62), An acoustic response signal corresponding to the acoustic wave is also 
measured (box 63). The measuring of the acoustic response signal (box 63) may 
be performed before the measuring of the seismo-electromagnetic signal (box 62), 
after the measuring of the seismo-electromagnetic signal (box 62), or 
simultaneously. 

[00142] The method further comprises exciting a formation with an electromagnetic 
exciting field (box 64). An electromagneto-seismic signal produced by the 
electromagnetic exciting field within the formation . is subsequently measured (box 
65). An electromagnetic response signal corresponding to the electromagnetic 
exciting' field is also measured (box 66). The measuring of the electromagnetic 
response signal (box 66) may be performed before the measuring of the 
electromagneto-seismic signal (box 65), after the measuring of the 
electromagneto-seismic signal (box 62), or simultaneously. 

[00143] The exciting with the acoustic wave (box 61), the measuring of the seismo- 
electromagnetic signal (box 62) and the measuring of the electromagnetic response 
signal (box 63) may be performed after the exciting with the electromagnetic 
exciting field (box 64) and subsequent measuring steps (boxes 65 and 66). More 
generally, the exciting steps (boxes 61 and 64) and the measuring steps (boxes 62, 
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63, 65 and 66) may be performed following any order that allows proper 
measurements. 

[00144] Once the seismo-electromagnetLc signal, the acoustic response signal, the 
electromagneto-seismic signal and the electromagnetic response signal are 
measured, an analysis may be performed so as to evaluate characterizing 
parameters of the formation. 

[00145] The acoustic response signal produced by the acoustic exciting (box 61) is 
an acoustic wave reflected at the formation or having traveled through the 
formation. Acoustic properties of the formation, e.g. velocities of various modes 
of acoustic waves, an acoustic impedance may be estimated from the acoustic 
response signal (box 67). 

[00146] The estimating of the acoustic properties (box 67) may for example be 
performed following an iterative inversion method: an initial value of the acoustic 
properties is selected. A synthesis acoustic response signal is synthesized from the 
initial value of the acoustic properties. An acoustic difference between the 
synthesis acoustic response signal and the measured acoustic response signal is 
calculated. The value of the acoustic properties is adjusted according to the 
acoustic difference. The synthesizing-, the calculating of the acoustic difference 
and the adjusting may be repeated until the acoustic difference drops below a 
predetermined acoustic threshold. 

[00147] Similarly, the electromagnetic response signal produced by the 
electromagnetic exciting (box 64) is sensitive to electromagnetic properties of the 
formation. Electromagnetic properties of the formation, e.g. an overall electrical 
conductivity a((o) or dielectric properties, may be estimated from the 
electromagnetic response signal (box 68). 

[00148] The estimating of the electromagnetic properties (box 68) may for example 
be performed following an iterative inversion method: an initial value of the 
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electroiiiagnetic properties is selected. A synthesis electromagnetic response signal 
is synthesized from the initial value of the electromagnetic properties. An 
electromagnetic difference between the synthesis electromagnetic response signal 
and the measured electromagnetic response signal is calculated. The value of the 
electromagnetic properties is adjusted according to the electromagnetic difference. 
The synfliesizing^ the calculating of the electromagnetic difference and the 
adjusting may be repeated until the electromagnetic difference drops below a 
predetemained electromagnetic threshold. 

[00149] Altematively, the estimating of the acoustic properties (box 67) and the 
estimating of the electromagnetic properties (box 68) may be performed with a 
method distinct from the iterative inversion method. 

[00150] Furthermore, the acoustic properties, e.g. a value of the velocities for 
longitudinal waves and shear waves, and/or the electromagnetic properties may be 
estimated from another logging measurement. 

[00151] The estimating of the acoustic properties (box 67) may also be performed 
after the estimating of the electromagnetic properties (box 68), or simultaneously. 

[00152] The estimating of the acoustic properties (box 67) and the estimating of the 
electromagnetic properties (box 68) allow to evaluate the characterizing 
parameters (box 69). The characterizing parameters may for example be an 
electrokinetic coupling coefficient X, a formation factor F of the formation and a 
mobility, i.e. a ratio k/ti of a permeability of the formation over a viscosity of a 
fluid within the formation. 

[00153] The permeability k may easily be deduced from the ratio k/ti and from the 
viscosity J/. 
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.[00154] The electrokinetic coupling coefficient L and the permeability k generally 
vary with a frequency of the acoustic wave or of the electromagnetic exciting 
field. 

[00155] FIG. 7 illustrates an example of an algorithm for evaluating characterizing 
parameters of a formation surrounding a borehole, according to an embodiment of 
the present invention. The algorithm of FIG. 7 may be executed by processing 
means of a system according to the present invention. 

[00156] Initial values of inversion parameters are selected (box 71). Typically, the 
inversion parameters are an electiroldnetic coupling coefficient Z and a mobility, 
i.e. a ratio k/ti of a permeability of the formation over a viscosity of a fluid within 
the formation. However, the inversiofn. parameters may be any other parameter 
involved in equations relating to electrokinetic phenomena. 

[00157] A synthesis seismo-electromagnetic signal and a synthesis electromagneto- : 
seismic signal are synthesized (box 72) using the initial values of the electrokinetic \ 
coupling coefficient Z, i.e. a coupling initial value, and the initial values of the 
ratio k/ii^ i.e. a mobility initial value. 

[00158] A first difference Dl is calculated between the synthesis seismp- 
electromagnetic signal and a measured seismo-electromagnetic signal and a 
second difference D2 is calculated between the synthesis electromagneto-seismic 
signal and a measured electromagneto-seismic signal (box 73). The calculating 
may for example involve a least squares method. 

[00159] The first difference Dl is compared to a fnst predetermined threshold Tl 
and the second difference D2 is compared to a second predetermined threshold T2 
(box 74). If the first difference Dl and the second difference D2 are respectively 
below the first predetermined threshold Tl and the second predetermined 
threshold T2, it is considered that the values of the inversion parameters are 
properly estimated. 
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[00160] If not, the values of the inversion parameters, i.e. the ratio k/n of the 
permeability of the formation over the viscosity and the value of the electrokinetic 
coupling coefficient L are adjusted according to the first difference and to the 
second difference (box 75), 

[00161] Hie synthesizing of the synthesis seismo-electromagnetic signal and of the 
synthesis electromagneto-seismic signal, the calculating of the first difference and 
of the second difference, and the adjusting are repeated until the first difference 
Dl and the second difference D2 drop respectively below the first predetermined 
threshold Tl and the second predetermined threshold T2. 

[00162] The synthesizing of the synthesis seismo-electromagnetic signal and of the 
synthesis electromagneto-seismic signal, the calculating of the first difference and 
of the second difference, and the adjusting may be performed following any order 
allowing a proper estimating of the inversion parameters. 

[00163] When the electrokinetic coupling coefficient L and the mobility, i.e. the 
ratio k/ti are considered as properly estimated, further characterizing parameters 
may be estimated. Typically, a value of a formation factor F of the formation may 
be determined (box 76). 

[00164] The determining of the formation factor may require to estimate the 
electrokinetic coupling coefficient L(a)) at different angular frequencies. A 
critical angular frequency a>c of the formation may be estimated asymptotically 
from the following equation, also established by Steven Pride (1994) : 

[00166] wherein : 

[00167] O represents an interconnected porosity; 
[00168] D represents an electrochemical potential; 
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[00169] Qoo represents a tortuosity; 

[00170] represents a pore fluid permeability; and wherein : 

[00171] a>, = (eq.4) 

[00172] When the critical angular frequency is known, as the ratio k/ti and the 
fluid density pf are also estimated, the formation factor F may be determined as : 

[00173] (eq.5) 
[00174] Le. F= • (eq. 6) 

[00175] The algorithm for evaluating the inversion parameters may be applied to the 
seismo-electromagnetic signal and to the electromagneto-seismic signal 
simultaneously, as represented in FIG. 7. The seismo-electromagnetic signal and 
the electromagneto-seismic signal are inverted simultaneously to evaluate the 
electrokinetic coupling coefficients and the mobility k/tj. 

[00176] Alternatively, one or more steps of the algorithm may be applied to a single 
signal and be subsequently repeated. 

[00177] The whole algorithm may even be applied to the electromagneto-seismic 
signal first, and then to Ihe seismo-electromagnetic signal, or vice-versa, before an 
additional step in which the electrokinetic coupling coefficient L and the mobility 
k/fi are estimated.firom intermediate parameters. 

[00178] Preferably the synthesizing takes into consideration a propagating of ati 
acoustic wave in the formation, in particular, the propagating of the acoustic wave 
in a matrix of a porous medium may be considered. The synthesizing is based on 
Pride's coupling equations (eq. 1 and 2). 
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[00179] Various modes of acoustic waves are generated within the formation, e.g. a 
fast longitudinal acoustic wave, a slow longitudinal acoustic wave, a shear 
acoustic wave, a Stoneley acoustic wave. Typically, in a case of an acoustic 
excitation : 

[00180] - the fast longitudinal acoustic wave produces a seismo-electric fast 
longitudinal signal only; 

[00181] - the slow longitudinal acoustic wave produces a seismo-electric slow 
longitudinal signal only; 

[00182] - the shear acoustic wave produces a seismo-magnetic shear signal 
only; 

[00183] - the Stoneley acoustic wave produces both a seismo-electric Stoneley 
signal and a seismo-magnetic Stoneley signal. 

[00184] The fast longitudinal acoustic wave and the slow longitudinal acoustic wave 
indeed fail to generate a magnetic field and the shear acoustic wave fails to 
generate an electric field. 

[00185] Similarly, in a case of an electric excitation, an electro-seismic fast 
longitudinal signal, an electro-seismic slow longitudinal signal and an electro- 
seismic Stoneley signal are produced within the formation. 

[00186] The synthesizing of a complete synthesis seismo-electromagnetic signal, i.e. 
a synthesis signal comprising a synthesis seismo-electromagnetic fast longitudinal 
signal, a synthesis seismo-electromagnetic slow longitudinal signal, a synthesis 
seismo-electromagnetic shear signal and a synthesis seismo-electromagnetic 
Stoneley signal, involves relatively complex calculation when based on Pride's 
coupling equation (eq, 1), 

[00187] The synthesizing may be simplified by considering only a single mode of 
acoustic wave, e.g. the slow longitudinal acoustic wave. 
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[00188] Alternatively, a plurality of well-chosen modes of acoustic waves is 
considered, e.g. the slow longitudinal acoustic wave and the fast longitudinal 
acoustic wave. 

[00189] Alternatively, a pluraKty of modes of acoustic waves may be separately 
considered, thus involving a plurality of relatively simplified calculations. 

[00190] Similarly, the synthesizing of a complete synthesis electromagneto-seismic 
signal involves relatively complex calculation when based on Pride's second 
coupling equation (eq. 2). The synthesizing may be simplified by considering only 
a single mode of electromagneto-seismic signal, e.g. the electro-seismic slow 
longitudinal signal. 

[00191] Alternatively, a plurality of well-chosen modes of electromagneto-seismic 
signals is considered, e.g. the electro-seismic slow longitudinal signal and the 
electro-seismic fast longitudinal signal. 

[00192] Alternatively a plurality of modes of electromagneto-seismic signals is 
considered separately, thus involving a plurality of smiplified calculations. 

[00193] The slow longitudinal acoustic wave has a relatively high sensitivity to the 
mobility. In tiie case of an acoustic excitation, the slow longitudinal acoustic wave 
hence is preferably considered. A synthesis seismo-electromagnetic slow 
longitudinal signal is synthesized from the mobility initial value and from the 
coupling initial value. 

[00194] The high sensitivity to the mobility of tiie slow longitudinal acoustic wave 
is well known from the art. The method of the present invention considers the 
seismo-electromagnetic slow longitiidinal signal associated to the slow 
longitudinal acoustic wave to provide information about the mobility. Unlike the 
slow longitudinal acoustic wave which is relatively difficult to measure directly. 



27 



PATENT APPLICATION 
ATrOKNEY DOCKET NO. 29244/030EPI 



the seismo-electromagnetic slow longitudinal signal may be measured relatively 
easily. 

[00195] Similarly, a synthesis electromagneto-seismic slow longitudinal signal is 
synthesized. 

[00196] The first difference is calculated between fte synthesis seismo- 
electromagnetic slow longitudinal signal and a slow longitudinal portion of the 
measured seismo-electromagnetic signal. 
[001971 The second difference is calculated between the synthesis electromagneto- 
seismic slow longitudinal signal and a slow longitudinal portion of the measured 
electromagneto-seismic signal. 
[00198] The first difference and the second difference are subsequently compared to 
predetermined thresholds, as described in FIG. 7. The values of the inversion 
parameters, i.e. the permeability of the formation over the viscosity and the value 
of the electrokinetic coupling coefficient are adjusted according to the first 
difference and to the second difference. The synthesizing steps, the calculating 
steps and the adjusting steps may be repeated, thus providing an iterative inversion 
method for characterizing the formation. 
[00199] The fest longitudinal wave, the Stoneley wave and the shear wave may also 

be considered, either alone or in combination with the slow longitudinal wave. 
[00200] Other simplifications may also be involved in the evaluating of the 

characterizing parameters. 
[00201] Alternatively, a distinct method may be used to analyze the seismo- 
electromagnetic signal and the electromagneto-seismic signal. 
[00202] The terms "seismic" and "acoustic" may be used indifferently in the present 
description: both terms relates to a mechanical phenomenon. 
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[00203] While the mvention has been described with respect to a Umited number of 
embodiments^ those skilled iu the art, having benefit of this disclosure, will 
appreciate that other embodiments can be devised which do not depart from the 
scope of the invention as disclosed herein. Accordingly, the scope of the 
invention should be limited only by the attached claims. 
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Claims 

A method for characterizing a formation, the method comprising: 

exciting the formation with an acoustic wave propagating into the formation (61); 

measuring a seismo-electromagnetic signal produced by the acoustic wave within 

the formation (62); 
exciting the formaition with an electromagnetic exciting field (64); 
measuring an electromagneto-seismic signal produced by the electromagnetic 

exciting field within the formation (65); 
analyzing the measured seismo-electromagnetic signal and the measured 

electromagneto-seisnuc signal to evaluate characterizing parameters of the 

formation (69)* 

The method of claim 1 , wherein 

the acoustic wave (408) and the electromagnetic exciting field (411) are generated 
at a logging tool (407; 507) positioned within a borehole (402; 502) 
sxirrounded by the formation (40 1 ; 50 1). 

The method of claim 1 or 2, further comprising: 

measuring an acoustic response signal, the acoustic response signal being 

produced by the acoustic exciting (63); 
estimating acoustic properties of the formation from the acoustic response signal 

(67); 

measuring an electromagnetic response signal, the electromagnetic response signal 

beinjg produced by the electromagnetic exciting (66); 
estimating electrofnagnetic properties of the formation from the electromagnetic 

response signal (68), 

The method of claim 3, further comprising: 
selecting initial values of inversion parameters (71); 
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synthesizing a synthesis seismo-electromagnetic signal and a synthesis 
electromagneto-seismic signal using the initial values of the inversion 
parameters (72); 

calculating a first difference between the synthesis seismo-electronaagnetic signal 

and the measured seismo-electromagnetic signal; 
calculating a second difference between the synthesis electromagneto-seismic 

signal and the measured electromagneto-seismic signal (73); 
adjusting the values of the inversion parameters according to the first difference 

and to the second difference (75); 
repeating the synthesizing using the adjusted values of the inversion parameters, 

the calculating of the first difference, the calculating of the second 

difference and the adjusting until the first difference and the second 

difference respectively drop below a first predetermined threshold and a 

second predetermined threshold. 

[cS] The method of claim 4, wherein : 

the inversion parameters are an electrokinetic coupling coefficient and a mobility; 

the synthesizing is simplified by synthesizing only a synthesis seismo- 
electromagnetic slow longitudinal signal and a synthesis electromagneto- 
seismic slow longitudinal signal from a mobility initial value and firom an 
electrokinetic coupling coefficient initial value. 

[c6\ The method according to any one of claims 1 to 5, wherein 

the analyzing takes into consideration the propagating of the acoustic wave (408) 
within the formation (401). 

[cT] The method according to any one of claims 1 to 6, wherein : 

the seismo-electromagnetic signal (409) is a seismo-electric signal. 

[cS] The method according to any one of claims 1 to 6, wherein : 
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the seismo-electromagnetic (409) signal is a seismo-magnetic signal. 

[c9] The method according to any one of claims 1 to 8, wherein 

the electromagneto-seismic signal (413) is a magneto-seisiiiic signal. 

(clOJ The method according to any one of claims 1 to 8, wherein 

the electromagneto-seismic signal (413) is an electro-seismic signal. 

[cll] The method according to any one of claims 2 to 1 0, further comprising: 

displacing the logging tool (407; 507) along the borehole (402; 502) so as to 
provide a continuous characterizing of the formation (401; 501) as a 
function of depth. 

[cl2J A system for characterizing a formation (401) surrounding a borehole (402), the 
system comprising: 

a logging tool (407) to be lowered into the borehole; 

an acoustic emitter (403) located onto the logging tool, the acoustic emitter 

allowing to excite the formation with an acoustic wave (408) propagating 

within the formation; 
an electromagnetic receiver (404) to measure a seismo-electromagnetic signal 

(409) produced by the acoustic wave within the formation; 
an electromagnetic emitter (404) located onto tiie logging tool, the electromagnetic 

emitter allowing to excite the formation with an electromagnetic exciting 

field (411); 

an acoustic receiver (403) to measure a electromagneto-seismic signal (413) 
produced by the electromagnetic exciting field within the formation; 

processing means (414) to analyze the measured seismo-electromagnetic signal 
and the measured electi-omagneto-seismic signal so as to evaluate 
characterizing parameters of the formation. 



32 



PATENT APPLICATION 
ATTORNEY DOCKET NO. 29244/030EP1 



[cl3] The system of claim 12, wherein : 

the electromagnetic receiver (404) is an electric receiver allowing to measure a 
seismo-electric signal (409) produced by the acoustic wave (408) wifliin the 
formation (401). 

[cl4] The system of claim 12, wherein : 

the electromagnetic receiver (504) is a magnetic receiver allowing to measure a 
seismo-magnetic signal produced by the acoustic wave within the 
formation (501). 

[cl51 The system of any one of claims 12 to 14, wherein : 

the electromagnetic emitter (404) is an electric emitter allowmg excite the 
formation (401) with an electric exciting field. 

[cl61 The system of any one of claims 12 to 14, wherein : 

the electromagnetic emitter (504) is a magnetic emitter allowing excite the 
formation (501) with a magnetic exciting field. 

[cl7] The system of claim 12, or 1 3 , further comprising: 

At least one additional electromagnetic receiver (506); 
At least one additional acoustic receiver (505). 
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Abstract 

CHARACTERIZING PROPERTIES OF A GEOLOGICAL 
FORMATION BY COUPLED ACOUSTIC AND 
ELECTROMAGNETIC MEASUREMENTS 

[00204] A method for characterizing a formation. The method comprises exciting 
the formation with an acoustic wave propagating into the formation. A seismo- 
electromagnetic signal produced by the acoustic wave within the formation is 
measured. The method further comprises exciting the formation with an 
electromagnetic exciting field. An electromagneto-seismic signal produced by the 
electromagnetic exciting field within the formation is measured. The measured 
seismo-electromagnetic signal and the measured electromagneto-seismic signal 
are analyzed to evaluate characterizing parameters of the formation. 
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Fig. 6 



Excite a formation with an acoustic wave propagating into the 
formation 



: IZ 

Measure a seismo-electromagnetic signal produced by the 
acoustic wave within the formation 
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Measure an acoustic response signal corresponding to the 
acoustic wave 
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Excite the formation with an electromagnetic exciting field 



IL 

Measure an electromagneto-seismic signal produced by the 
electromagnetic exciting field within the formation 

, V 

Measure an electromagnetic response signal corresponding to the 
electromagnetic exciting field 

_iz — — 

Estimate acoustic properties of tiie formation firom the acoustic 
response signal 

iz 

Estimate electromagnetic properties of the formation fix>m the 
electromagnetic response signal . 

\/ ..^ 

Estimate characterizing properties of the formation firom the 
seismo-electromagnetic signal and firom tite electromagneto- 
seismic signal. 
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